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ABSTRACT 

A wind tunnel simulation of the diffusion patterns in a sea breeze 
has been attempted. No attempt was made to reproduce the recircu- 
lation that characterizes a sea breeze, but the results indicate 
that the low level onshore flow was well simulated for neutral, 
stable, unstable, and elevated inversion conditions. Velocity, 
turbulence, shear stress, and temperature data were taken, and 
the spread of emissions from ground level sources was investigated. 

Comparison is made with theoretical predictions by E. Inoue and 
with the open, homogeneous plane field results of Pasquill. 
Agreement with the predictions by Inoue is good. The comparison 
with Pasquill’s results shows that the wind tunnel flows are 
shifted two categories towards more stable. The discrepancy may 
be explained as a matter of averaging time. 



2 




J. INTRODUCTION 


Near large bodies of water, sea breezes ( or lake breezes) are 
common features of the climatology. The different surface temperature 
and roughness between sea and land causes complicated stratification 
of the air flow, which tends to create an elevated inversion. This 
causes fumigation and plume trapping (Lyons and Cole, 1973). This 
and the occurrence of return flow make, the sea breeze an important 
factor when studying the dispersion of pollutants near lake or sea 
shores where large industries are commonly situated. 


Figure 1 is a schematic diagram showing the flow patterns commonly 
occurring in a sea breeze. The height of the boundary layer above 
the sea may be 300m to AOOm The sea breeze flows from 5km to 
as much as 40km inland, and then rises, forming a front. The 
flow then reverses and returns seaward at a height of 600 m or 
more, usually not exactly 180° to the onshore flow.’ It may then 
gradually descend and flow onshore again, making a complete cycle. 
Within the sea breeze, the velocity and temperature profiles 
change continuously since the roughness and temperature on the 
land surface are different from those of the sea surface. 


Yoshida, Aoki and Takamatsu (1972) found that the sea breeze advances 
as a wedge near the ground and the diffusivity may become very small 
in this wedge. Fukuoka (1968) investigated high concentrations 
occurring upwind from a stack. Such occurrences may be explained by 
the return flow which feeds pollutants to the sea breeze as it invade 
the land. Lyons (1970) observed a similar return of pollutant at lake 
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Michigan. Wakamatsu, et al. (1971) found S0 2 concentrations behind 
a sea breeze front five to ten times higher than the concentrations 
before the sea breeze began. After the sea breeze died, concentrations 
returned to the morning values. Nakano (1972) observed a similar vari- 
ation of concentration during the sea breeze cycle. 

Vith current techniques it is not yet possible to simulate the complete 
sea breeze, including return flow in the wind tunnel. The experiment 
reported here is an attempt to simulate the surface layer flow and to 
Study the gross flow characteristics — mean velocity, turbulence, and 
temperature profiles — and the concentration dispersion patterns for 
ground level sources located on the land under the influence of ther- 


mal stratification. 



2. THE WIND TUNNEL MODEL 


k 


The typical surface boundary conditions of the sea breeze might 
be described as follows: 


sea: small roughness length (z QS ), lower surface temperature (0 S ) 

land: large roughness length (z 0i j), and/or elevated topography, 
higher daytime surface temperature (Q^). 

These conditions also describe suburban and urban areas respectively, 

,so similar conditions may be used to simulate both sea breeze and urban 

heat island effects on diffusion. 


The turbulent boundary layers in the actual and model sea breezes may 
be compared using the turbulent Reynolds number. and densimetric Froude 
number. The turbulent Reynolds number at height H is 

Re tU rb= J (H) . H , 

kTh) 


where K is the turbulent diffusion coefficient (eddy di ffus ivi ty) . If 
a logarithmic wind profile is assumed, then K(H) = ku^H and 


Re 


turb ~ 


= i ln ( |L) 

ku^H k 2 ° 


which does not depend on the wind velocity but only on the roughness 


length. 


Generally, z Q and the boundary layer thickness, 6, over the sea are about 
10" 1 cm and 350 m respectively; therefore, Re tlirb - 80. If 6 = 15 cm in 
the wind tunnel, then z Q * 5 x 10~ 5 cm is required to simulate the above 
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In a flow field with both mechanical and thermal turbulence, the 
two Influence each other. A surface to air temperature gradient 
changes turbulent components and these Influence the mean velocity 
profile and turbulent diffusion coefficient. These in turn alter 
the temperature profile. Thus, the temperature and velocity fields 
are not independent. In this sense, Richardson number, a ratio of 
mechanical to thermal energy, becomes important for diffusion pheno- 
mena in thermally stratified fields. 


Four surface temperature conditions were considered in this study, 
corresponding to neutral, unstable, stable and elevated inversion 
conditions. A maximum of A0 = 30° C was chosen as the sea- land 
surface temperature difference. This is small compared with the 
required temperature difference from Froude number, however, a large 
vertical temperature gradient above the ground created Richardson 
numbers the same order of magnitude as might be expected in the 
atmosphere, although for short downwind distances. 


The gradient Richardson number is given as 

30 

Ri =*' g 3z 

Mi ) 2 

Minus values indicate unstable conditions, plus values indicate 
stable conditions. If J RI J is nearly zero, mechanical turbulence 
dominates, whereas a large [Ri | shows that thermal turbulence dominates 


Qf POOR QUALTOS 
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The relation between L (Monin - Obukov length) and Ri can be shown 
to be: 


Ri 


z a 

i- 


where a is the ratio of turbulent diffusion coeffients of momentum 
and heat and 0m is a nondimen sional shear term. 


The mean velocity profile for near neutral conditions is given by: 



which is the so called log-plus-linear law, where $ is a constant. 

For unstable conditions, L is less than zero, making plots of log z 
vs U concave upwards. For stable conditions, L is greater than zero, 
so the semi -log velocity curves become concave down. 

If the sea is colder than the air above it, the air over the water 
becomes stably stratified. As this air moves over the land, it is 
heated from below by the land, and the stratification near the ground 
changes to unstable. Therefore ,the vertical profile of the Richardson 
number changes from minus to plus with increasing height. The height 
where Ri = 0 corresponds to the height of the base of the inversion 
which develops. Thus the Richardson number profile can be important 
and useful for describing the flow and diffusion fields. 
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3.„ DIFFUSION UNDER NON-NEUTRAL CONDITIONS 
Many of the presently used diffusion equations do not appropriately 
describe diffusion in the atmosphere under non-neutral conditions. 
Empirical theories, such as Pasquill’s, work well under non-neutral 
conditions over a homogeneous plane, but are difficult, if not 
impossible, to apply in the case of a sea breeze. Such theories, 
designed for constant conditions in x and y, typically cannot account 
for the longitudinal variation of stability that occurs in a sea breeze. 


The results obtained from this experiment will be compared with several 
useful relations developed by E. Inoue (i960). I noue considered three 
conditions: extremely unstable, neutral and extremely stable and found 

the following approximate relations for o y , and the ground level con* 
centration as functions of downwind distance: 


extremely unstable a 


neutral 


z«Q 

°Y 

a z 

C z =0 


ec 'X 
* 


1,5 


.- 2*5 


0-8 


r— 1*8 


extremely stable 




z=0 


0.5 


,-1 • 5 


These relations may be used in a step-by-step method in comparison with 
experimental results for flows with changing stability. As stability 
changes with distance from the shoreline, the appropriate power relation 
changes. Log-log graphs of the experimental results can therefore be used 
to indicate the stability category and the approximate strength of the stability 
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7 . k. EXPERIMENTAL APPARATUS AND MEASUREMENT TECHNIQUE 
The tests were performed in the Air Pollution Wind Tunnel at New 
York University. The wind tunnel test section is J.l m x 2.2 m x 10 m 
(3-5 ft. x 7 ft. x 35 ft.) and has heaters for controlling the air in- 
take temperature and maintaining the same wall, floor and ceiling 
temperatures. A survey carriage allows the remote positioning of 
probes to the nearest 0.25 mm. 

The model for this experiment consisted of six smooth aluminum plates 
providing 365 cm of “sea" and 365 cm of "land". The sea portion was 
raised 2.5 cm from the wind tunnel floor so that the step would gener- 
ate a turbulent boundary layer). The land portion was raised an addi- 
tional 1.25 cm for topographic similarity, as discussed earlier. Twelve 
heating tapes with separate controllers were fastened to the underside 
of the land surface to provide temperatures up to 100° C with good 
uniformity. The sea surface was cooled by slipping slabs of dry ice 
Into the space between the aluminum plates and the wind tunnel floor. 
Flush mounted, 2.5 cm diameter sources covered with stainless steel mesh 
were added to the land at downwind distances of 2.5* 25 and 100 cm from 
the Step between the sea and land. Sixteen thermocouples were fastened 
to the upper surface to be used for recording and adjusting the surface 
temperature. The resulting configuration is shown in Figure 2 . 

* 

Velocity data were measured and recorded using a DISA constant tempera- 
ture anemometer system. Mean velocity, vertical and longitudinal turbu- 
lent velocities, Reynolds stress and heat flux were determined by a pro- 
cedure described by Arya (1968). This technique involves rotating a 



single hot wire , so a probe rotater was constructed to turn the probe 
support about its long axis in precise 45° Intervals. A fine gage 
thermocouple and a resistance thermometer were also fastened to the 
rotater to measure temperature profiles and the r.m.s. of the temper- 
ature fluctuations. 

The sources built into the land were fed a mixture of 2 % ethane in 
nitrogen. This hydrocarbon tracer was sampled from a single brass 
probe on the survey carriage using a Beckman hydrocarbon analyzer. 

“Hie output was recorded and time-averaged with one minute averaging. 
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5. DISCUSSION OF RESULTS 
5“1* Flow Measurements ' 

The following graphs summarize the flow data in nondimensionalized 
form for the four test conditions with crossplots for comparison. 

All distances are nondimensionalized by the step height (d 51 1.25 cm) ^ 
x* =* x/d is the nondimensional downwind distance from the step at the 
sea-land interface, and z* = z/d is the local nondimensional height 
above the plate% except in the plot of log z* vs. U* (Figure 3d) f where 
z* is measured from sea level. Wind speed, turbulent velocities and 
Reynolds stress are nondimensionalized using the free stream velocity 
Uos = 1.0 m /s. Temperatures are plotted in Centigrade degrees over or 
under wind tunnel ambient, 30° C. The labels Neutral, Unstable, Stable 
and Inversion are applied in a generic sense to simplify the description 
of the wind tunnel conditions. 

The flow characteristics were measured for all four stability categories, 
but complete data for only the Inversion case are presented here for the 
sake of brevity. Figure 3 shows the nondimensionalized mean velocity, 
the mean temperature and the Richardson number profiles for the Inver- 
sion case. The results Illustrate the characteristics of a typical ele- 
vated inversion. The velocity profiles in Figure 3a show the rapid 
readjustment which occurs downwind of a step in a flow field. In Figure 
3b the temperature gradient is unstable below the inversion base and 


stable above the base. 
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The Richardson number profiles show the inversion details more clearly. 
Richardson number was calculated using the data contained in Figures 
3a and 3b and the gradient form 

30 

R1 * g 3z 

w" 

The Richard numbers have a wide range, so a logarithmic scale has been 
used in Figure 3c with an overlap at Ri ** ± 10"^. The height where 
Ri =* 0 corresponds to the elevation of the inversion base. The height 
of the inversion base increases almost linearly with downwind distance. 
The semi-log velocity profiles. Figure 3d, have the characteristics pre- 
dicted by the log-plus-linear law of Monin and Obukhov. In the unstable 
layer below the inversion, the profiles appear concave down and reverse 
to be concave upward in the stable layer above the inversion base. 

Using the inflection points where the concave up and down portions meet 
as an indication of the inversion height, it may be seen that this predic- 
tion yields approximately the same heights as does the Richardson 
number prof* le. 

Figure 4 shows comparisons of the Neutral, Unstable, Stable and Inversion 
flow data at x* - 20. In Figure 4a, which shows mean velocity profiles, 
it may be seen that the Neutral and Unstable boundary layer thicknesses 
are approximately the same, but the Unstable profile shows increased 
velocity at lower height due to increased vertical turbulence, which 
augments momentum transfer from higher to lower levels. The Stable and 
Inversion boundary layer thicknesses are smaller, with less uniform 
velocity. 
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Figure Ab shows the corresponding mean temperature profiles for 
the Unstable, Stable and Inversion conditions. The unstable layer 
thickness is z* < The stable layer thickness is z* < 8, consid- 
erably greater than the unstable layer because the stable layer was 
developed over an additional 3&5 071 of upwind r, sea". The Inversion 
case temperature profile is unstable for z* < 2 and stable for 
2 < z* < 8. 

Longitudinal turbulence ( ) , vertical turbulence ( yw 2 ”) , and 

Reynolds stress (-uw) are smallest for the Stable case flow and greatest 
for the Unstable case, as shown in Figures k c, d, and e. Neutral cases 
values typically fall between the Stable and Unstable values. Inversion 
case values are close to the Stable for 2 < z* < 8, and closer to the 
Unstable case values for z* < 2. Figure kf shows that the vertical 

heat flux is a small minus value for the Stable case flow and a large 
plus value for the Unstable- flow* • The heat flux for the Inversion case 
is a small minus value in the stable layer above the inversion base and 
a large plus value in the unstable layer below the base. 
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5“2. Concentration Measurements 

Concentration measurements were carried out for the four air stabili- 
ties (Neutral, Unstable, Stable and inversion) and three different 
ground level source locations. For each of these twelve cases, ground 
level crosswind profiles and vertical profiles were determined at five 
downwind distances from the source. Distances have been nond i mens ional i zed 
by the step height as per the velocity data. The sources are designated 
simply No* 1, No. 2, and No. 3 at nond i mens iona 1 distances of 2, 20 and 
80, respectively from the step (shoreline). Except where the omission 
would neglect important details only data concerning Source No. 2 is 
presented here. For each source, signifies the distance from the 
center of the source, in order to simplify comparison. The concentra- 
tions were normalized by the concentration measured at the downwind 
edge of- the 2.5 cm. diameter source. These are multiplied by 10 3 for 
convenience. 

Figure 5 shows the vertical concentration profiles at four downwind posi- 
tions from Source No. 2 for the four stability conditions. In general, 
the Unstable profiles have the greatest spread ,fo1 lowed by Neutral, 

Inversion and Stable. On the other hand ground level concentrations gen- 
erally were highest for Stable, followed by Inversion, Neutral and Unstable 
as a consequence of the different spreads. 

Figure 6 shows the vertical profile at x* ; - 20 for Source No. 2 with 
greater detail. The Neutral case profile has almost a linear decrease 
with increasing height for the semi-log coordinates shown indicating an 
exponential decay. The Unstable case profile also appears exponential. 
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- however, another point off the graph ( the trend is indicated by the 
line ) and the profiles at other positions show that the curve is 
not exponential. The Stable and Inversion case profiles are both 
considerably compressed; the Inversion case is slightly less so 
because of the unstable layer below the inversion base. These 
compressions correspond to the decreased vertical turbulence near 
the ground as shown in Figure 

Vertical standard deviations (a^) , defined as the second moment of 

the profile about its mean were calculated for each profile. The 

values of . - o z were nond Imens ionalized by the step height and are 

presented in Figure 7» Data from Source No. 1 have been included 

because they show details which became obscured downwind of Source 

No. 2, where the inversion became weaker. The values of are 

largest for the Unstable condition, and decrease for Neutral and 

Stable. The three slopes indicated at the bottoms of the graphs 

correspond to the predicted increases for extremely unstable ( x 3/2 ) , 

neutral (x 0 * 8 ) and stable ( x 1/2 ) „ The experimental Stable results 

have a smaller slope than predicted, and the Unstable and Neutral 

have slightly smaller slopes than predicted. Since Source No. 1 

was near the step, where heating had just begun, it may be expected 

that the unstable stratification would be weak, and the stable 

would be strong, because the source was close to the sea. The new 

internal boundary layer created by the step might be expected to 

affect o^, but the changes due to stability appear to' be more 

significant. It is interesting to note that the Inversion a has 

z 

almost the,, same value as the Stable , but Increases rapidly to 
the Neutral and Unstable values. Since the flow is predominant ly 
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stable near Source No. 1, and the unstable layer underneath becomes 
thicker with distance, such a rapid increase is reasonable. 

Figure 7b shows that the characteristics for Source No. 2 were 
generally the same as for Source No. 1, with the expected changes 
that accompany the greater distance. The trapping of the plume by 
the inversion layer can be seen more clearly in Figure 8, which 
shows the heights at which the concentrations equal one-tenth of 
the peak concentrations for each source under the inversion condition. 
The heavy dashed line corresponds to the height of inversion base 
according to the Richardson number profiles. The plume boundary 
height from Source No. 1 is almost the same as the inversion height 
which suggest trapping due to the stable layer above this height. 

The plume boundary for the Neutral case is higher close to the source 
but lower farther downwind. The plume boundaries from Source No. 2 
are almost the same under Neutral and Inversion conditions, while 
the height under Inversion conditions is greater than under neutral 
for Source No. 3. At these downwind locations, the unstable stratifi- 
cation under the inversion raises the plume boundary higher than it was 
for the neutral case. 

Figure 9 shows ground level crosswind profiles at four downwind 
distances from Source No. 2 under the Neutral, Unstable, Stable 
and Inversion conditions. As already shown by the vertical profiles 
the magnitude of the ground level concentrations on the centerline 
depends on stability with the same relationship of Stable as highest 
and Unstable lowest. 
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The spreads also follow the same pattern: widest for Unstable 
and Neutral and narrowest for Stable. 


The crpsswind profiles can be accurately described by normal 
distributions. The standard deviations a were determined by 

y 

plotting the nond imensional concentrations on probability paper, 
from which the standard deviation and location of the mean can be 
read easily. Because 2,5 cm diameter sources were used in this 
study, does not fall to zero at x** « 0, but had a value of 
tfy s 0.75 cm. This Is reasonable if a normal distribution is 
accounted at the source position. To convert the area source data 
to point source data for comparison with field studies, “ 0.75 cm. 
has been subtracted from the wind tunnel values. These values of 
0 ^ have then been nondimens ionalized using the step height to obtain 
a y* “ °Y /d - 


Figure 10 shows plots of the log of these converted <jy* values versus 

x*, for Sources No. 1 and 2. In each case, the curves show that 

a y^ * xl > approximately. Unstable values of are largest, followed 

by Neutral and Stable, the same relations as for a z *. For Source No. 

the Inversion Oy* also begins at a small value, the same as Stable, 

and increases to a value greater than the Neutral value., However, 

for Source No. 2, the Inversion are much closer to the Unstable 

values, as was the case for a *. 

z 


Figure 11 shows the centerline ground level decays of concentration 
with downwind distance for Sources No. 1 and 2. According to Inoue 
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(1969), the ground level concentrations may be proportional to 
x “ 1 » 5 "for extremely stabl e, >x _1 • 8 for neutral, and x for 

extremely unstable. These slopes are drawn in the figure for 
comparison. For both sources, the agreement between predicted and 
experimental results is good except for Source No. 1 at the smaller 
values of x*' , which may be due to step effects. 
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5*3. Comparison with Field Diffusion Experiments 

Most recent field experiments have been compared with Pasquill's six 
stability categories, and the agreement for the available experimental 
data over flat, homogeneous planes is quite good. But near the sea- 
shore, where sea or land breezes occur frequently, these curves will 
not always fit the results of field diffusion experiments. 

Using data from their field experiment near the seashore, Sakuraba 
et al. (1967) plotted calculated values of over Pasquill's curves. 
These were shifted two categories in the stable direct ion, and , if the 
longer sampling time were accounted for (30 min. vs. 10 min. for Pas- 
quill's curves), the shift would become even more significant. 

Other field studies conducted near the sea show a gradual shift of 
a z from stable category F to neutral category D with distance. Such 
results show that near the seashore the flow is not homogeneous in the 
longitudinal direction as it is over open, flat land, where Pasquill's 
curves were established. 

To compare this wind tunnel data with Pasquill's curves, scaling is 
required. Scales of 1000:1, 2000:1 and 3000:1 were investigated when 
recalculating the wind tunnel data for comparison with Pasquill's curves. 
At 1000:1, the data points spread from x = 125 to 1000 meters. Increas- 
ing the scale ratio increased this distance and shifted stability class 
from stable towards neutral. The effect, however, was small, so 2000:1 
will be considered an acceptable, representative scale for this comparison 

Figure 12 shows a z vs. x for source No. 3 at the scale ratio 2000:1 for 
the four stability classes tested. The data extends over the range 250 
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to 2000 meters. The Stable condition curve falls into category F, 

Neutral into E and Unstable Into C to 0. At the source location, 

Inversion values lie between Neutral and Unstable due to the unstable 
surface layer. 

Figure 13 shows cr y vs. x for Source No. 3 at the scale ratio 2000:1 
for the four stability classes. Wind Tunnel Neutral fits F, Unstable 
fits E, and Stable is much more stable than F. The wind tunnel limits 
horizontal eddy size, which increases with observation time in the 
field. Because of the plume meandering the mean concentrations decrease 
with increased sampling time. 

Hino (1966) proposed a 1/2 power decay of with observation time. 

This was determined from turbulent diffusion theory and agrees well 
with experiment for long sampling time (several minutes to several hours). 

The wind tunnel data, when multiplied by 2.5» agree well with Pasquill's 
curves. As Pasquill's curves are for 10 minutes samples^the 1/2 power 
law then implies that the equivalent wind tunnel sampling time is 1.6 
minutes or 96 seconds. This means that the wind tunnel sampling time 
corresponds to a small observation time for diffusion phenomena in the 
field . 

The normalized centerline ground level concentration used by Pasquill is: 

C U = i 
-G * a y a z 

For ground level releases and measurements like those made In this study 
the value of U is difficult to define. Since a y and a z were measured in 
the wind tunnel, the normalized concentrations were estimated from these. 
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Figure 1 A shows the normalized ground level concentration, based on 

the calculated a and a , plotted over Pasquill's curves. A scale 
Y 2 

ratio of 2000:1 and Source No. 3 were used for the comparison. The 
results again show a shift of two categories towards stable. 
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6. CONCLUSIONS AND FUTURE STUDY 

Sea breeze effects on diffusion were simulated in the wind tunnel 
by heating and cooling the wind tunnel floor. Flow and diffusion 
characteristics were measured for four stabilities. Close to the 
shoreline, stable or inversion conditions increased the ground 
level concentrations behind the ground level sources to more than 
five times the neutral values. Since most wind tunnel diffusion 
studies have been carried out for neutral stability, this result 
Is especially significant. Inflow of stable air from a sea or 
lake is a common part of the climatology for seaside industrial 
parks, so stability effects due to sea breezes are important to 

air quality. 


In the present study, diffusion from a ground level source was 
considered. For a more realistic study, an elevated source model 
Is desirable. Such a study will be next in the series of thermal 
experiments being conducted by the authors. The relation of source 
height to inversion layer height would be very important in a 
study of fumigation, where higher ground level concentrations may 
be expected than would occur for an elevated source with neutral 
conditions. 
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Sketch of commonly occuring sea breeze flow patterns. 

Schematic of the model. The vertical scale is exaggerated 
to show the step. A summary of the boundary temperature 
conditions is given beneath the sketch. 

Summary of the downwind development of the measured flow 
characteristics for the inversion condition. 

(a) Nondimensional ized mean velocity profiles 

(b) Mean temperature profiles 

(c) Richardson number profiles, shewing inversion height 
(dashed 1 ine) 

(d) Semi-log mean velocity profiles 

Comparison of the measured flow characteristics for each 
stability condition. 

(a) Mean velocity profiles 

(b) Mean temperature profiles 

(c) Longitudinal turbulent intensity 

(d) Vertical turbulent intensity 

(e) Nondimens ional ized Reynolds stress 

(f) Vertical heat flux 

Summary of the downwind development of the vertical 
concentration profiles for each stability condition. 

Source No. 2. 

Comparison of the vertical concentration profiles at 
x* /= 20 for the four stability conditions. 

Source No. 2 

Vertical standard deviations of the plume spread vs. 
downwind distance from the source. 

(a) Source No. 1 

(b) Source No. 2 

Plume width for the Neutral and Inversion conditions vs. 
downwind distance. 

Summary of the development of the ground level crosswind 
concentration profiles for the four stability conditions. 
Source No. 2. 

Crosswind standard deviations of the plume spread vs. 
downwind distance from the source. 

(a) Source No. 1 

(b) Source No. 2 
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Figure IK Decay of the center line ground level concentration 
with downwind distance from the source. 

(a) Source No. 1 

(b) Source No. 2 


Figure 12. Comparison of wind tunnel results for cr with Pasquill’s 
curves. The wind tunnel data have been mu! t ipl ied by a 
scale ratio of 2000:1, as explained in the text. 


Figure 13. 


Comparison of wind tunnel results for a with Pasquill’s 
curves. The wind tunnel data have been^mul t i pi ied by a 
scale ratio of 2000:1 , as explained in the text. 


Figure 14. Comparison of normalized ground level concentration in 
the wind tunnel with Pasquill’s curves. The wind tunnel 
data have been multiplied by a scale ratio of 2000:1, 
as explained in the text. 
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4 f) VERTICAL HEAT FLUX 




NEUTRAL 



INVERSION 



Figure 6 




Figure 7a 


n » 


3.0 


7 b) SOURCE NO. 2 


O NEUTRAL 



Figure 7b 



O PLUME EDGE , NEUTRAL 
A PLUME EDGE , INVERSION 













Figure 10b 
















